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ABSTRACT 

I investigate whether useful constraints on the evolution of the relationship 
between galaxy mass (jrigai) and black hole (BH) mass (ttibh) can be obtained 
from recent measurements of galaxy stellar mass functions and QSO bolometric 
luminosity functions at high redshift. I assume a simple power- law relationship 
between mgai and mBH, as implied by BH mass measurements at low redshift, 
and consider only evolution in the zero-point of the relation. I argue that one can 
obtain a lower limit on the zero-point evolution by assuming that every galaxy 
hosts a BH, shining at its Eddington rate. One can obtain an upper limit by 
requiring that the number of massive BH at high redshift does not exceed that 
observed locally. I find that, under these assumptions, and neglecting scatter in 
the mgai-mBH relation, BH must have been a factor of ~ 2 larger at 2; ~ 1 and 5-6 
times more massive relative to their host galaxies at z ~ 2. However, accounting 
for intrinsic scatter in TOgai-mBH considerably relaxes these constraints. With 
a logarithmic scatter of 0.3-0.5 dex in mBH at fixed mgai, similar to estimates 
of the intrinsic scatter in the observed relation today, there are enough massive 
BH to produce the observed population of luminous QSOs at 2; ~ 2 even in 
the absence of any zero-point evolution. Adopting more realistic estimates for 
the fraction of galaxies that host active BH and the Eddington ratios of the 
associated quasars, I find that the zero-point of the mgai-mBH relation at 2; ~ 2 
cannot be much more than a factor of two times larger than the present-day 
value, as the number of luminous quasars predicted would exceed the observed 
population. 
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1 INTRODUCTION 

The discovery of the relationship between the masses of 
nuclear supermassive black holes (SMBH) and the lumi- 
nosity, veloci ty disper si on, or m ass of their host galactic 



spher oids ( Dresslej 



I995I: Magorrian et al 



198 9|: i Kormendv fc Richstone 



200d: iGebhardt et al.1 



1199. 



iFerrarese fc Merritt 



2OOOI : iMarconi fc Hunt! l2003l : 



Haring fc Rixl 120041 ) is surely one of the most profound 



observational results of the past decade, if not the past 
century. Different methods applied to both dormant 
and active BH in the nearby Universe now yield con- 
sistent results and indicate that BH mass and galaxy 
velocity dis persion a are related via mBH oc , where 
p ~ 4-5 (IFerrarese fc MerrlttI I2OO0I : IGebhardt erahl 
I2OO0I : iTremaine et al.ll2002l ). and BH mass and galaxy 
mass are related via rriBH oc rrigai arconi fc HuntI 
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l2003l : iHaring fc Rb3l2004l fl The observed scatter in both 
of these relationships is remarkably small, and implies 
intrinsic scatters of approximately 0.3 de x in ttibh at 
fixed g, and 0.5 dex in t hbh at fixed L (|Novak et al.l 
l2006l : iGultekin et alJlioogl ). 

A large number of theoretical explanations for 
the origin of this observed relationship (hereafter re- 
ferred to for brevity as the mBH-ro-gai re l ationship) have 



been proposed (e.g . ISilk fc Reed Il998l: iBurkert fc Silkl 
1 lAdams etaiT l200ll. l2003l: IWvithe fc Loebl 1200^ 



2001 



Robertson et al.1 l2006l : ICrotonI l2006l : [Hopkins et al.l 



200 7a'). However, there is no widely accepted unique the- 
oretical model, and the models differ in their predictions 
for the amount of evolution in the mBH-T7^gal relationship. 
This quantity is therefore a potentially strong discrimi- 
nator between different theoretical models, and there is 
great interest in obtaining robust direct observational 
measurements of this relationship at high redshifts. A 
great deal of effort and telescope time has been expended 
towards achieving this goal. 

With currently available facilities, the masses of dor- 
mant BH can be measured via the dynamics of th e sur- 
rounding gas or stars (see iFerrarese fc FordI [20051 . for a 
recent review) only for very nearby galaxies. At high 
redshift, it is currently possible to attempt to measure 
masses only for active or accreting BH. Several such stud- 
ies have claimed to find evidence for significant evolution 
in the mBH-'Tigai relationship, always in the sense that 
black holes are more massive at high redshift relative to 
their ho st galaxy (e.g. Treu et al. 2004; Pong et al. 2006"; 
IWoo et al.,,2006., 2008. ; Salviandor ct al. 2007). However, 
these methods rely on a set of simplified underlying as- 
sumptions and various proxies for the desired quanti- 
ties, and are sub j ect to potentially severe selection bi- 
ases. ILauer et al.1 (|2007al ) have argued that if there is a 
moderate amount of scatter in the intrinsic mBH-Wgai re- 
lationship (consistent with observational constraints on 
the scatter in the local relation), these selection biases 
can account for most or all of the claimed evolution. It is 
therefore interesting to explore the possibility of obtain- 
ing independent empirical constraints on the evolution. 
Even obtaining robust upper and lower limits on the evo- 
lution could be useful. 

Recently, it has become common practice to estimate 
the stellar masses of galaxies from multi-wavelength 

broadband photorn e try and/ or sp ectroscopy (e.g. 

Bell fc de Jon3 12OOII : [Bell et all l2003l : iKauffmann et all 



20031 ). The availability of deep multi-wavelength surveys 



covering substantial sky area has yielded estimates of the 
stellar mass f unctio n s out to high r e dshift , z ~ 4-5 (e.g. 



Drorvetal.1 |2004|: iFontana et all |2004|: lBorch etal.1 



20061: 



Font ana et al.l 



2006: iPanneUa et al.l 



-r 



20061: 



Bundv et al.' '2006'; 'Pozzctti ct al.' '20o7; "Vergan i et al, 
2008; iMarchosini ct al. .2008; P erez- Gonzalez et al, 
20081 ). Similarly, there has been significant progress 



^ Several other equally tight relationships between BH mass 
and other galaxy properties have been disco vered (e.g. 
iGraham et~^ I2OOII : iKormendv fc Bended |2009| ). Although 
these relationships are intriguing, here I focus on the rela- 
tionship between BH mass and galaxy mass. 



in piecing together the evolution of the bolometric 
luminosity function of quasars and Active Galactic 
Nuclei (AGN) out t o high redshift (z '-^ 6) from multi- 
wavelength surveys (iHopk ins et al ] |2007bl . and references 
therein; hereafter HRH07). 

In this paper, I explore whether one can derive useful 
empirical constraints on the evolution of the relationship 
between BH mass and galaxy mass by comparing these 
two sets of observed statistical distributions (galaxy stel- 
lar mass functions and QSO/AGN luminosity functions) 
under the basic and well-accepted ansatz that accreting 
BH provide the power source for AGN. The outline of the 
rest of the paper is as follows. In <}21 I describe my ba- 
sic set of assumptions. In i|2.1I I present upper and lower 
limits on the evolution of the mBH-m-gai relation, assum- 
ing that there is no intrinsic scatter in the relation. In 
i]2.2I I discuss how these constraints are impacted by the 
inclusion of intrinsic scatter. I conclude in 33] 



2 UPPER AND LOWER LIMITS ON THE 
EVOLUTION 

For nearby dormant BH, the average relationship be- 
tween BH mass a nd galaxy mass can be characterized 
as m BH oc TTiga/'^ (|Marconi fc Huntj[200^ : iHaring fc Rfa3 
[200313. I will explore the simplest possible form for pos- 
sible evolution of this relationship, namely scaling by a 
factor that is a function of redshift only r{z). Thus the 
mass of a BH residing in a galaxy with mass rrigai at 
redshift z is given by: 

mBH(2;, rrigai) = r{z) mBH{z = 0, mgai) = r(2:) m^^i^'^ . 

(1) 

I now make use of the observed galaxy stellar mass 
function at some redshift of interest, and assume that 
every galaxy hosts a SMBH with mass given by Eqn. [1] 
One can then obtain a reasonable lower limit on the evo- 
lution in mBH-Wgai (i.e. on ^{z)) at a given redshift by 
assuming that 1) every BH is active at all times (has a 
duty cycle of unity) 2) every active BH always radiates 
at its Eddington luminosity. This set of assumptions will 
clearly maximize the number of luminous quasars for a 
given population of BH, under the fairly standard con- 
jecture of Eddington-limited accretion. 

One can then obtain an upper limit on the evolu- 
tion by comparing the implied BH mass function at the 
redshift of interest with observational estimates of the 
present-day BH mass function. Again, under the appar- 
ently reasonable assumptions that BH masses increase 
monotonically with time and that significant numbers of 
massive BH are not somehow lost from galaxies, clearly 
the number of massive BH in the past cannot exceed that 
at the present day. To state the condition more precisely, 
the number density of BH, (j){mBu), at high redshift may 
not exceed the present day value for all BH masses greater 
than some threshold value man > Afmin. 

For the observed stellar mass functio n as a function 
of redshift, I adopt the fitting functions of iFontana et al.l 



^ iLauer et al.l ll2007bl) suggest a slope of unity, but this would 
not significantly efi^ect the results. 
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Figure 1. The bolometric luminosity function of quasars at ^ = 1 (left panel) and z = 2 (right panel). Square symbols with 
error bars show the estimate of the observed bolometric QSO LF from HRH07. Dashed lines show the upper limit on the QSO 
LF derived from the observed stellar mass function at the relevant redshift and the arguments described in the text. From left 
to right, the dashed lines assume that the zero-point of the mBH-"*gal relation has evolved by a factor of F = 1, 2, 3, 4, 5, or 6. 
Under these assumptions, the mass of a typical BH hosted by a galaxy of a given mass must have been larger by a factor of ~ 2 
at 2 = 1 and by a factor of 5-6 at 2 = 2. 



l|2006l ). based on measurements from the GOODS- 
MUSIC survey. I have checked that these fitting func- 
tions provide good ag reemen t with the results from other 
surveys (see Fontan ot et al.l l)2009h for a thorough com- 
parison of stellar mass function estimates from different 
surveys over a wide range of redshifts). There is good 
agreement between different estimates of the stellar mass 
function up to 2: ~ 2; at higher redshifts the results of 
different studies diverge. For this reason, the results pre- 
sented here are limited to redshift two and below. The 
stellar masses have been converted to correspond to a 
IChabriej l|2003l ) stellar initial mass function. I will assume 
that there are no significant redshift-dependent system- 
atic errors in the stellar mass estimates. For the moment, 
because I am mainly exploring the feasibility of this ap- 
proach, I also ignore the effect of the random errors on 
the stellar masses, although these will probably have an 
impact on the quantitative results. 

2.1 Constraints without Scatter 

Initially, I assume that the relationship between BH 
mass and galaxy mass has no intrinsic scatter. It is then 
straightforward to derive the implied bolometric lumi- 
nosity function of QSO/AGN from the observed galaxy 
stellar mass function under the set of assumptions out- 
lined above, for a given value of the evolution factor T{z). 
Figure [T] shows the observed bolometric QSO luminosity 
function as estimated by HRH07 along with the upper 
limit estimate based on the stellar mass function, for dif- 
ferent values of P. This comparison is shown ai z — 1 and 
z = 2. At QSO luminosities below the "knee" in the LF, 
the upper limit estimate overproduces QSOs, which can 



be understood as implying that these BH are not active at 
all times and/or are radiating at sub-Eddington luminosi- 
ties. What is more interesting is that at the highest lumi- 
nosities, above ~ at z = 1 or 10^'^'^ at z = 2, with- 
out evolution in the mBH-'Tigai relation, the number of 
luminous QSOs is significantly underestimated. Assum- 
ing that these luminous QSOs are not radiating above 
their Eddington luminosity, and are not magnified some- 
how (e.g. by beaming or lensing), one can then read off 
the minimum amount of evolution (minimum value of F) 
required to produce enough luminous QSOs. This corre- 
sponds to Fmin ~ 2 at 2: = 1 and Fmin ~ 5-6 at z — 2. 
Taken at face value, then, these results require that BH 
were 5-6 times larger for a given galaxy mass at z = 2. 

Now let us consider the upper limit on the evolution, 
or maximum allowed value of F. Figure [2] shows the ob- 
servational estimate of the BH mass functioijfl at z = 
(|Marconi et al.|[2004l ). compared with the results from the 
stellar mass function at 2: = 2 scaled by the same series 
of values of F. Clearly, as long as BH cannot decrease in 
mass or be ejected from their host galaxies, then there is 
an upper limit of Fmax ~ 6 at z ~ 2. It is interesting that 
the lower limit from the QSO LF, discussed above, and 
this upper limit are so close to one another. 



3 Note that the lMarconi et al.l ||2004) estimate of the BHMF 
includes the effect of scatter in the local mBH-™.gal relation, 
and therefore the comparison shown here is not strictly self- 
consistent. However, I will consider scatter self-consistently in 
the next Section. 
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Figu re 2. The mas s funct ion of SMBH. Tlie solid green line shows the observational estimate of the BH mass function at z = 
from lMarconi et alj l l2004h . Dashed (red) lines show the BH MF implied by the observed galaxy stellar mass function at z = 1 
(left panel) or z = 2 (right panel), the relationship between BH mass and galaxy mass described in the text, and evolution in the 
zero-point of the "igH-^gal relation of a factor of F = 1, 2, 3, 4, 5, or 6 (curves from left to right, respectively). 




Figure 3. Left: BH mass function. Open (purple) dots show the BHMF implied by the observed z = 2 GSMF, no evolution in 
the m,BH-"^gal relation (F = 1), and a scatter of (tbh = 0.3; err or bars are simple Po isson errors. The solid (green) line shows 
the observational estimate of the BH mass function at z = from lMarconi et al. I 1I2OO4I) : the dashed (red) line shows the BH MF 
implied by the observed galaxy stellar mass function at z = 2 plus the assumed mBH-"^gal relation with no evolution (F = 1) and 
no scatter. To satisfy the constraint, the purple dots (prediction) should be lower than the green line (observations). Right: The 
QSO luminosity function. Open (purple) dots show the upper limit on the QSO LF from the same argument, but including scatter 
in the »nBH-™gal relation (ctbh = 0.3); error bars are Poisson. Square (green) symbols with error bars show the estimate of the 
observed bolometric QSO LF at z = 2 from HRH07; the dashed (red) line shows the upper limit on the QSO LF derived from the 
observed GSMF at z = 2 and the arguments described in the text, for F = 1 and ctbh = 0. Here, to satisfy the constraint, the 
purple dots (prediction) should be higher tha.n the green squares (observations). The inclusion of a moderate amount of scatter 
has a large impact on the high-mass end of the BHMF and the high luminosity end of the QSO LF. When scatter is included, the 
assumption that the rraBH-"*gal relation has not evolved since z ~ 2 appears to be consistent with these constraints. 



© 0000 RAS, MNRAS 000, 000-000 



Evolution of the mBH-?^gai Relationship 5 




o 



6 7 8 9 10 10 11 12 13 14 15 

log [Mq] log L^o, [L©] 



Figure 4. The same as Fig.fS] but with ctbh = 0.5. 
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Figure 5. The same as Fig. |3] but with F = 2 and (tbh = 0.3. 



2.2 Constraints with Scatter 

So far I have neglected intrinsi c scatter in the mBH -w2gai 
relation. However, eis noted bv iLauer" et al.l l|2007ah . be- 
cause of the very steep decline of the high-mass end of 
the galaxy mass or luminosity function, the highest mass 
black holes are actually more likely to be outliers from 
the mBH-?Tigai relation, hosted by galaxies of more modest 
mass, rather than typical BH in the much rarer high-mass 
galaxies that would be the sole hosts of such BH in the 
absence of scatter. 



iNovak et al.l (|2006| ) attempted to constrain the in- 
trinsic scatter in the mBH — cr and mBH — L relations at 
z — 0, and concluded that due to the small sample of 
galaxies with reliable measurements of BH mass and un- 
certainties in the observational error estimates on ttibh, 
only upper limits on the scatter could be obtained. They 
estimated these upper limits to be Sa < 0.3 and 5l < 0.5, 
where Sa is the 1-a log scatter in the mB H — o" relation and 
5l IS t he same for the ttibh — L relation. Im arconi fc HuntI 
(|2003l ) find a simil ar scatter in the mBH — msph relation 
as m man — o", and lHaring fc Rixl (|2004l ) find an observed 
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scatter of 0.3 dex in the rriBH — Jrigai relation, implying 
that the intr insic scatter is presumably smaller. Recently, 
iGultekin et al . (2009) made a detailed study of the mag- 
nitude of the intrinsic scatter in mBH — o" and ttibh — L, 
finding 5o- = 0.31 for ellipticals but a larger scatter of 
5a — 0.44 for all galaxies (including spirals). They fur- 
thermore find that the shape of the distribution of the 
intrinsic residuals in rriBH at fixed a is consistent with a 
log-normal (and inconsistent with a normal distribution) . 

Unfortunately, almost nothing is known about the 
possible evolution of the intrinsic scatter in the BH- 
galaxy scaling relations. Therefore, I investigate how the 
inclusion of representative amounts of scatter would im- 
pact the results presented in Section [2.11 In order to do 
this, I run Monte Carlo simulations of ~ 10^ galaxies, in 
which I first select galaxy masses from the observed stel- 
lar mass function at the redshift of interest (using z = 2 
as a representative case) . I then assign BH masses to each 
galaxy according to Eqn. [1] adding a random deviate in 
mass selected from a log-normal distribution with root 
variance ctbh, and consider the implied BH mass func- 
tion and upper limit on the QSO LP as before. 

Results of these experiments for various values of 
(tbh and F (all at z = 2) are shown in Figures |3H5] 
Note that I cut off the galaxy stellar mass function be- 
low lO^^M© because these low-mass galaxies do not pro- 
vide interesting constraints and including them causes 
the Monte Carlo simulations to take much longer to run 
(for a given desired number of high-mass objects). In Fig- 
ure O one can see that when a moderate scatter in the 
niBH-m^ai relation is included {aBH ~ 0.3, similar to the 
scatter in the observed relation at the present day), the 
number of luminous QSOs can be reproduced under the 
fairly extreme assumptions used in the lower limit exer- 
cise (all BH radiate at their Eddington limit at all times) . 
Even a scatter half as large as the observed present-day 
estimates {aBH ~ 0.15), with no evolution in the normal- 
ization (F = 1) marginally satisfies the lower limit. With 
a slightly larger scatter {aBH ~ 0.5) or moderate evolu- 
tion in the zero-point F = 2 (see Fig. |4]and[5]), bright 
QSOs are overproduced by a factor of 10-100, leaving 
room for more reasonable assumptions about duty cycle 
and Eddington ratio. 

One can try to sharpen this constraint by adopt- 
ing more physically reasonable va lues for the dut y cy- 
cles and Ed dington ratios of AGN. lErb et al.l (|2006D and 
iKriek et all (2007) find that about 20-40% of galaxies at 
2: ~ 2 contain an active nucleus, and models in which 
such activity is merge r-driven (e.g. iHopkins et al. I l2008l : 
ISomerville et al.1 l2008i) predict that this fraction is nearly 
constant for galaxy masses 10.0 ^ log(m « /Mo) ^ 12.0 

iseee.g. Figure 19 of iHopkins et al" I I2OO8I ). IVestergaardI 
20041 ) find that the Eddington ratios of luminous quasars 
at 1.5 < 2 < 3.5 are in the range 0.1 < L/L-Edd < 
1, with an average value L/L^dd ~ 0.4-0.5, while 
iKoUmeier et al.l ([2006') find a fairly sharply peaked distri- 
bution of Eddington ratios with a peak at L/L^dd ~ 0.25. 
Adopting average values for the fraction of galaxies con- 
taining active BH, /agn ~ 0.3, and the Eddington ratio 
/Bdd = L/LEdd — 0.5 (assuming that L/L^dd is also con- 




10 11 12 13 14 15 
log Lbol [L©] 

Figure 6. The QSO LP at 2 = 2 as shown in Fig.[3l with F = 1 
and (Tbh = 0.3, but with a QSO duty cycle of /agn = 0.3 and 
an Eddington ratio of /ecM = 0.5. 

slant with galaxy/BH massjf] produces quite good agree- 
ment with the observed QSO LF at 2 = 2 with no evolu- 
tion in the zero-point or scatter of the mgai-mBH relation 
(F = 1, asH = 0.3; see Fig.©. 

Turning the argument around, then, if the indepen- 
dent observational estimates of duty cycle and Eddington 
ratio are correct, and if the scatter in the mgai-mBH rela- 
tion was not significantly smaller at high redshift than it 
is today, then overall evolution in the mgai-rriBH relation 
of F 2 since z ~ 2 is disfavored as it would overproduce 
the number of luminous QSOs (see Fig.[71l. In particular, 
if the value of F at « ~ 2 were as large a s suggested by the 
observations of e.g. I Peng et al.l (|2006l ). F ^ 4, the num- 
ber of luminous QSOs would be overproduced by more 
than one order of magnitude. Of course, one could recon- 
cile these larger amounts of evolution if the duty cycle of 
luminous quasars is an order of magnitude smaller than 
what I have assumed (~ 2-3 percent instead of 20-30 
percent). 



3 CONCLUSIONS 

I have investigated whether observational estimates of the 
stellar mass function of galaxies, combined with observed 
QSO luminosity functions, can provide useful limits on 
the relationship between galaxies and their SMBH at high 
redshift. I assumed a simple relationship between galaxy 
mass and SMBH mass, as observed in dormant galaxies 
in the nearby Universe, and a simple form for the possible 

* Note that if we included a realistic distribution of Eddington 
ratios, rather than a single constant value, this would again 
broaden the tail of the bright end of the QSO LF, leading to 
more very luminous QSOs. 
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evolution of this relationship (see Eqn.[T|, namely a shift 
in the zero-point of the relation by a redshift-dependent 
factor r(2). I then argued that one can obtain a lower 
limit on V{z) by making the rather extreme assumption 
that all BH radiate at their Eddington limit at all times, 
and requiring that at least the observed number of lu- 
minous QSOs be reproduced. I further argued that an 
upper limit on T{z) could be obtained by requiring that 
the number of massive BH in galaxies today should not 
be exceeded at high redshift. 

Assuming that there is a deterministic relationship 
between galaxy mass and BH mass (i.e., no scatter in the 
'TiBH-'Tigai relationship), I find that in order to produce 
enough luminous QSOs, the zero-point of the relation 
must have been higher by at least a factor of ~ 2 at 2 = 1 
and a factor of 5-6 at z = 2. At the same time, in order to 
avoid producing a larger number density of massive BH 
than what is implied by observations at 2 ~ 0, the upper 
limit on the evolution of the normalization of the ttibh- 
mgai relationship at 2 = 2 is about a factor of six. Since 
both the lower and upper limits are fairly liberal, one 
might have expected them to lie several orders of magni- 
tude apart, and therefore not to provide very interesting 
constraints on the actual evolution of the mBH-Tigai rela- 
tionship. It seems potentially quite interesting that these 
limits lie nearly on top of one another. 

However, relaxing the assumption of a perfectly de- 
terministic mBH-Tigai relationship has a major impact 
on the results. When scatter is included in the rriBH-wigai 
relation at a level similar to the intrinsic scatter in the 
observed relation at 2 = 0, 1 find that the majority of very 
massive BH are objects that live in galaxies of moderate 
mass but are outliers in the mBH-?Wgai relationship. This 
is of course due to the very steep slope of the galaxy stel- 
lar mass function at large masses. Because the constraints 
above arose from the most luminous QSOs, I then find 
that there is a strong degeneracy between the evolution 
of the zero-point r(2) and the scatter (Jbh- For example, 
the QSO constraint at 2 = 2 can be reproduced even 
in a scenario in which F = 1 (no evolution in the zero- 
point has occured) and (tbh ~ 0.3 (the intrinsic scatter 
in mBH-?Tigai is similar to that in the observed relation 
today). Thus we are left with the very weak constraint 
that BH probably were no smaller at high redshift rela- 
tive to their host galaxies (unless the scatter was much 
larger than it is today). 

I tried to sharpen this constraint by adopting more 
physically reasonable values for the duty cycles and Ed- 
dington ratios of AGN, based on independent observa- 
tional constraints. Adopting mass-independent values of 
,/agn = 0.3 (the fraction of galaxies hosting AGN) and 
/Edd = L/L-Edd ~ 0.5, and assuming a scatter in niBH- 
mgai similar to that in the observed relation for dormant 
galaxies today (asH = 0.3), I find that BH cannot have 
been much more than a factor of ^ 2 more massive rel- 
ative to their host galaxies at 2 ~ 2 than they are to- 
day. In particular, values as large as r(2 — 2) ~ 4, as 
sug gested by some observational studies (e.g. IPeng et al.l 
I2OO61 ). would overproduce the number of luminous QSOs 
by more than an order of magni t ude. 

Interestinglv. lHopkins et al.l (|2006h also reached sim- 
ilar conclusions based on a somewhat different, though re- 



lated argument. They pointed out that in order to avoid 
overproducing the total mass density in SMBH relative 
to the present day value, the average value of mBn/jrigai 
must not have been more than about a factor of two larger 
at 2 ~ 2 than today's value. 

I have based these results on the relationship be- 
tween the total stellar mass of the galaxy and the mass 
of the SMBH; however, there is strong evidence that 
the more fundamental relationship is actually between 
the BH mass and the mass of the spheroidal compo - 
nent of the galaxy (e.g. iKormendv fc Richstond Il995l ). 
I made this choice because the stellar mass function of 
galactic spheroids is very poorly constrained at high red- 
shift. However, at low redshift, the most massi ve galaxies 
are p redominantly spheroid-dominated (e.g. iBell et ahl 
I2OO3I ). If this was also the case at high redshift, then 
my conclusions will not change much as the constraints 
are driven by the most massive BH which are hosted by 
massive galaxies. If there is a significant population of 
disk-dominated massive galaxies at high redshift, and the 
BH mass indeed correlates with spheroid mass only, then 
this would leave more room for evolution and/or scatter 
in the rngph-m-BH relation. 

Another source of uncertainty arises from the fact 
that BH masses predicted from the ttibh vs. luminos- 
ity {mBH-L) relationship are inconsistent with those pre- 
dicted from the ttibh vs. velocity dispersion (msH-o") re- 
lations hip for the most luminous galaxies l|Lauer et al.l 
l2007bl ). The rriBH-JTigai relationship that I have chosen 
to use here is deriv ed from the tubh-L relation, which 
iLauer et al.1 (l2007bl ) argue should be more reliable in the 
regime of interest, but the situation at high redshift is 
unknown. Currently, there are no published observational 
measurements of the galaxy velocity dispersion function 
at high redshift (of which I am aware); however, these 
may become available in the future. It would then be very 
interesting to repeat this kind of analysis using mBH-o" 
instead. 

Although it is dissappointing that the proposed ap- 
proach did not yield stronger constraints on the evo- 
lution of the mBH-'Tigai relationship, this exercise has 
brought out a few important lessons. First, in order to 
understand the relationship between galaxies and their 
BH, it is perhaps as important to understand the mag- 
nitude and evolution of the intrinsic scatter in this re- 
lationship as it is to understand the evolution of the 
zero-point of the relation itself. Second, new generations 
of theoretical models that attempt to simultaneously 
treat the formati on and evolution o f galaxies and their 



black holes (e.g. ICroton et al.|[2006l: iBower et all I2OO6I : 



iFontanot et al.n2006l : ISomerviUe et al.l I2OO8I ) must take 
care to properly model the dispersion in the mBH-?Tigai 
relationship. 
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Figure 7. The QSO LF at 2 = 2 as shown in Fig. |3] with intrinsic scatter ctbh = 0.3, a QSO duty cycle of /agn = 0.3 and an 
Eddington ratio of /fidd = 0.5. Left panel: F = 2; Right panel: F = 5. Assuming that the duty cycles and Eddington ratios derived 
from independent observations are correct, and that the intrinsic scatter in the rvigai-mBH relation was at least as large at z = 2 
as it is today, large amounts of evolution in the zero-point (F ^ 2) are disfavored. 
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